Excitations of the One Dimensional Bose-Einstein Condensates in a Random Potential 
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We examine bosons hopping on a one-dimensional lattice in the presence of a random potential 
at zero temperature. Bogoliubov excitations of the Bose-Einstein condensate formed under such 
conditions are localized, with the localization length diverging at low frequency as l{u) ~ l/u a . 
We show that the well known result a = 2 applies only for sufficiently weak random potential. 
As the random potential is increased beyond a certain strength, a starts decreasing. At a critical 
strength of the potential, when the system of bosons is at the transition from a superfluid to an 
insulator, a = 1. This result is relevant for understanding the behavior of the atomic Bose-Einstein 
condensates in the presence of random potential, and of the disordered Josephson junction arrays. 
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One of the most challenging problems of quantum 
many body physics is the behavior of stongly interact- 
ing matter in a disordered environment. In this paper 
we investigate the universal properties of superfluids in 
such systems, near the superfluid insulator transition. In- 
terest in this problem arises in many independent con- 
texts, in work on granular superconducting films and 
wires [TJ [51 [3] , Helium condensates in vycor [3] , and re- 
cent experiments on Bose condensates in optical traps. 
In particular, issues such as the expansion of a noninter- 
acting Bose condensate through a random potential 
excitations in an interacting Bose Einstein condensate in 
a random potential [6j [7] , and the possibility of the ob- 
servation of the Bose glass phase [HI [5] were explored in 
very recent theoretical and experimental papers. As im- 
portant is the possibility of investigating the behavior of 
disordered superconductors in a controlled fashion using 
Josephson junction arrays, as in Refs. (TUJ [TTJ [12]. In 
low dimensional quantum systems, where symmetry bro- 
ken phases are very fragile, we expect the most dramatic 
manifestations of the interplay of disorder and interac- 
tions. The existence of the Bose-glass phase was estab- 
lished in Refs. [13j [14], where the scaling and renormal- 
ization group (RG) picture of the Id superfluid-insulator 
transition at weak disorder was also established. Re- 
cently, much theoretical progress was afforded through 
real-space RG approaches in the case of dissipative [TS] 
and closed [16j [17] bosonic chains, where the properties 
of the SF-insulator transition at strong disorder were es- 
tablished. 

In this paper, we study the excitations of the super- 
fluid phase in a bosonic chain with a strongly random 
potential and interactions, near the SF-insulator transi- 
tion. Capitalizing on the real-space RG understanding 
of this transition [TBI HZji we analyze the localization 
length of phonons (i.e., Bogoliubov quasiparticles) as a 
function of their frequency and wave number. Deep in 
the superfluid phase, when the random potential is weak, 
the phonon localization length t{ui) at small uj diverges 



as mrriinn 



a = 2. 



(1) 

(2) 



This result, in particular, formed the basis of the anal- 
ysis in Refs. [BJ [7] . Using the renormalization group anal- 
ysis of Ref. [TS] and the study of random elastic chains of 
Ref. [19], we show that Eq. ([2| does not apply everywhere 
in the superfluid regime. In a finite region of parameter 
space on the superfluid side near the superfluid-insulator 
transition, Eq. |2| fails, and is replaced by the law 



a 



(3) 



where 1 < g < 2. The meaning of the parameter g 
will be elucidated later in the paper. Furthermore, as 
the system approaches the transition to the insulating 
regime, g decreases. Exactly at the transition g = 1, and 
Eq. acquires a correction to scaling: 



t(uS) ~ (ln 2 w) juj. 



(4) 



Eqs. ([3| and Q constitute the main result of our paper. 

Our analysis begins by considering a one-dimensional 
disordered Bose-Hubbard model with many particles per 
site. Its Hamiltonian is 



Uk 
2 



d 



Jk cos ( 



k) 



(5) 



This Hamiltonian describes a chain of sites, connected to 
their nearest neighbors by a Josephson hopping with a 
random strength, Jfe. Uk is the strength of the onsite re- 
pulsion, and rife S [— 5,5] are random offset charges. 
The hopping, charging and offsets are randomly dis- 
tributed with probability densities Pj(J), Pu(U), and 
P n (n). 

In the strong-disorder limit a real space renormaliza- 
tion group analysis can be employed to gradually elimi- 
nate sites with anomalously large Jk or or local charging 



2 



gap, Afe = [7fc(l — 2|nfc|) [TBI EZ] ■ The remaining sites are 
described by the same Hamiltonian but with the renor- 
malized probability distributions. The system of Eq. ([5| 
then emerges as either a superfluid or an insulator; the 
latter could be either a Mott insulator, a Mott glass, Bose 
glass, or random-singlet glass, depending on the strength, 
relative and absolute, of various types of disorder present. 
If the bosonic system is a superfluid, the distribution of 
J renormalizes towards the universal limiting function 



Pj(J) = cj g - 



(6) 



with C providing normalization. The superfluid is de- 
scribed by g > 1 with its value decreasing as the critical 
point at g = 1 is approached; in particular, as disorder 
increases, g decreases. At the same time, 



Pu{U) 



IP 



exp 



u 



(7) 



where / flows to along the renormalization group tra- 
jectories and f2 is the decreasing UV cutoff scale of the 
renormalized Hamiltonian, i.e., its largest hopping or 
gap. We now proceed to show that the same parameter g 
appearing in the distribution ^ controls the localization 
length of low-frequency phonons, as expressed in Eq. Q3J) . 

At the final stages of the renormalization, as long as 
g > 1 , the system is a superfluid and the possibility that 
the phase difference at adjacent sites slips through 2ir 
can be safely ignored. Moreover, since the Hamiltonian 
is no longer periodic in <pk it is now possible to do a gauge 
transformation to eliminate rik in Eq. dSJ, marking the 
ability of the superfluid to screen arbitrary offset charges. 
Then one may expand the cosine, to find the effective 
Hamiltonian (up to an unimportant additive overall con- 
stant) 



Ukd^ 

2 ddn 



Jk 
Y 



-4>uY 
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The Hamiltonian Eq. ([8| is quadratic, and thus we obtain 
full information by analyzing it at the classical level. Its 
classical equations of motion are 



to 2 u: 



l 4>k = Jk - 0fe+i) + Jk-i {4>k - <t>k-x) , (9) 



where lo is the angular frequency. These describe 
phonons in a chain of random masses connected by ran- 
dom springs. The masses are m& ~ while the 
spring constants are proportional to Jk- Ref. [19] pre- 
sented the solution to this problem (referred to as "Dyson 
type II") for the case when P V (U) = 5{U - U ) (that 
is, nonrandom uniform charging energies Uk = Uq) and 
Pj(J) = CJ 9 ~ l with J € [0, J ]. It was found that the 
average density of states is constant at low frequency, 



p{u>) = ( 6 (uj — u> n ) J ~ const. 



(10) 



Here ui n are the frequencies of the phonons described by 
^ and brackets denote averaging over random Jk- A 
constant density of states at low frequency is of course 
a feature shared with non-random elastic chains. The 
phonons which are the solutions of Eq. ^ , however, are 
all localized. Their localization length obeys Eq. (JTJ (at 
small oj) with 



a = 2, 5 >2; a = g, 1 < g < 2. 



(11) 



These results are compatible with our claim in beginning 
of this paper, Eqs. ^ and Q. Nevertheless, the uniform- 
U treatment leading to Eqs. ( 11 1 cannot be considered a 



derivation of the localization length results in our prob- 
lem: the random boson problem has charging energies 
Uk which are also randomly distributed. Below we show 



that the results given in Eqs. ( 10 1 and ( 11 ) are valid even 



if Uk are random, as long as the probability of observing 
anomalously small Uk is not too large. In addition we 
show that the uniform- U results for g = 1 exhibit strong 



corrections to scaling [see Eq. (31 1]. 

Let us first confirm that the fully random Bosonic 
chain, Eq. also has a finite constant density of states 
at low frequencies, as in (10 1. Consider the classical 
ground state of a system of N + 1 sites described by 
Eq. pj, where the first grain's phase 0o is kept fixed at 
4>a = 0, and a force h is applied conjugate to the phase 
4>n oi the last grain. The equilibrium values of the vari- 
ables (f>k can be found by minimizing the energy 



JV-l 



E = | ^2 J k (4>k - 4*k+if ~ h(j) N . 



(12) 



k=0 



which yields: 



k-l 



(13) 



Alternatively, 4>k can be computed in the following way. 
Introducing the variables ipk = 4>k/VUk, we can find ipN 
by inverting the matrix Tiki defined by the expression 

I N ^ 2 i N 

2 £ Jk (^kVu~k- iPk+ix/Uk+ij = - ^2 i~ikiipki>i- 



k=0 



Then 



k,l=l 



(14) 



4>k = y/U k U N G kN h. (15) 
where G is a matrix inverse to TL. In particular, we are 



interested in k = N case when Eq. (15) can be rewritten 
as 



dn) 
'N 



N 

E 

fc=i 



u k 



(16) 
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(n) 

Here 4> k are the normalized solutions to the cigcnmodc 
equation Eq. ^ with the boundary conditions 4>o = at 
the beginning of the chain, and with the frequency u> n , 
labelled by the index n. 

Next we compare the two expressions for cf>k at k = N, 



Eq. ( 13 ) and Eq. ( 16 1. We observe that for the probabil- 
ity distribution P(J) = J 9_1 , as long as g > 1, 



IN-l 



>N 



N. 



(17) 



, fc=0 



On other other hand, 



4>n = N dui p(oj) 



a 



(18) 



where 4> refers to the eigenmode at frequency u>, and 
p(oj) is the density of states. Clearly, unless the density of 



states is strongly suppressed at small uj, the integral ( 18 ) 
diverges due to small ui contributions. At small u>, the 
localization length exceeds the system size, thus (fiffl ~ 
1/y/N. At the same time, (1/Uf.) is finite, which means 
that C u is both uj and N independent. Suppose p(ui) ~ 
w 7 , where 7 < 1. Then 



du> 



p(oj) p(uj ) 



(19) 



where uj is the smallest frequency of the system, which 
can be found by 



u « 1 1 

dujp(uj) ~ — : , UJ - 

M iV !+"> 



This in turn gives 



>N 



(20) 



(21) 



Comparison with (17) reveals that 7 = 0, i.e., Eq. < \10\ . 

We now return to the localization length. First, con- 
sider the case of weak random = J + SJ k and 
Uk = Ua + SUk- Treating SJk as a perturbation, it is easy 
to find the localization length following Refs. [5U1 |2"T] . 
Indeed, the mean free time can be found by the Fermi 



golden rule, to go as r 



while the mean free path 



goes as \J JoUqt. The localization length is proportional 
to the scattering length in ID, thus l(oS) ~ ui~ 2 . 

This calculation however ignores the possibility of the 
wave scattering off the anomalously small Jfc or Uk- In- 
deed, suppose we have a "weak link" in Eq. ([9} where 
J W eak link = j on that link is much smaller than J on 
other links, j -C J- It is straightforward to check that 
the phonons with wave vector q ^> j/J get reflected off 
this weak link, while those with wave vector q <C j/J 
pass straight through. This is easiest to see if we solve 
Eq. ([9]) with the assumptions that all Jk are equal to Jo, 
while that of the weak link is j -c Jo an d all the Uk are 



equal to Uq. Then the transmission coefficient through 
the weak link is given by 



T = 



1 



J 2 ' 

4j2 



(22) 



where q is the dimensionless wave vector which is as- 
sumed to be small, or \q\ <C 7r. T tends to 1 at small q, 
and to at large q » j/ J. 

It thus follows that a phonon with wave vector q can- 
not have a localization length bigger than the average 
distance between the "weak links" with the strength of 
their couplings no bigger than q divided by the density 
of states. Using Eq. ^ we can estimate the average sep- 
aration between such weak links. We find 



dJ J 9 ~ 



1 



(23) 



where £ is the average distance between the weak links 
j. This gives £ ~ l/j 9 - Since j ~ q, and q ~ uj due to 
Eq. (10} , 

1/cj 



by 



the localization length is bounded from above 
9 , thus we arrive at our result, Eq. 
Scattering off the small Uk can also reflect the short 
wavelength waves. Taking all Uk equal to Uq, while the 
"heavy link" (recall that the "masses" are inversely pro- 
portional to Uk) Uk equal to u <C Uq, and taking all the 
Jfc = Jo gives the transmission coefficient 



T = 



1 



(24) 



4u 2 ' 



equivalent to (22 1. This, however, does not lead to any 
corrections to Eq. ( 3|. Indeed, using the same arguments 
as preceeding Eq. ( 23 1 , we find 



dU 



Of 

u 



exp 



u J 



(25) 



Here £ is the typical distance between these "heavy" 
links. Again taking u ~ q ~ u>, we find 



exp 



(26) 



This estimate is much bigger than Eq. ^ and thus the 
real localization length Eq. ([3| remains unaffected. This 
concludes the derivation of Eqs. (JTo|) and Eqs. ( 2|3|4 1. 

The analysis of £(oj) above assumed that we probe the 
phonon modes of the superfluid in the very end of the 
renormalization group flow, once the power law that con- 
trols the distribution Pj(J) ~ J 9-1 , has already attained 
its fixed-line value. While this is valid in the limit of 
uj — > 0, corrections to scaling may arise near the critical 
point. Refs. [HI [17] allows us to consider the corrections 
to this analysis arising from the flow to the SF fixed line. 
The RG flow for the generic-disorder case is given by: 



df ,, , dg 1 



(27) 
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where T is the logartihmic RG flow parameter. In the 
region close to the critical point, / = Q, g = 1, we can 
solve these equations approximately to give: 



9 



1 



2c 



4e e 



(28) 



for disorder realizations that flow to g = 1 + e with e <C 1 . 
Flows that terminate at the critical point, however, are 
given approximately by: 



1 



1 f ~ A 



(29) 



To find the corrections to scaling in the form of £(uj), 
we first note that it is given by the bare length-scale of 
rcnormalized sites once the RG scale reaches T = In — , 
with uiq the bare energy scale of the Bose-Hubbard chain. 
The RG flow of the effective site and bond length is: 



dl 



(30) 



At the critical point we expect £(u)) ~ 1/w; let us first 
derive the correction to scaling at the critical point. In- 
tegrating Eq. ( |30| using Eq. (29 1 gives ln£ = 

2inr/r + 0(l/T), 



r 



using Eq. 

and thus we find the localization 
length at criticality having a logarithmic correction: 



£{u) ~ [In 2 (w/wq)] /w. 
Off criticality, we find by the same analysis: 

%)~[(i-( w /^)V*] 2 /w 1+e . 



(31) 



(32) 



In summary, localization properties at low frequency 
are determined by the parameter g. Its value cannot be 
calculated directly in closed form from the initial disorder 
distribution, but we can estimate it by following the RG 
flow using the techniques of Refs. [THl HI]. In Fig. [I] 
we demonstrate how initial distributions evolve into the 
exponent a, which is g at the end of the flow. 

Finally, we remark that the results of this paper should 
also be valid in the case of the quasi one-dimensional 
condensates in the presence of random potential (but in 
the absence of any lattice). Indeed, such condensates 
are expected to form puddles in the minima of the po- 
tential, with nonzero boson hopping amplitudes between 
the puddles. Then they are expected to mimic ([5]), and 
the rest of the analysis of this paper applies. 
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FIG. 1: An example of the correspondence between initial 
distributions, and the exponent g = a at the end of the 
flow. Main plot: terminal g as a function of the param- 
eter 8 for initial Gaussian coupling distributions Pu(U) ~ 
e -(ff-o.s)»/o.7» ) and Pj(J) „ e -(J-s) 2 /0A 2 t truncat ed for 
J,U < 10~ 4 . Inset: example of the flow of g vs. the RG 
flow parameter V for the Gaussian initial conditions with 
5 — 1.125. The localization exponent a at 5 — 1.25 is g at 
the large P plateau. Each point is averaged over 40 disorder 
realizations, with the chains of 5 • 10 6 sites. 
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